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Supercritical state has been viewed as an intermediate state between gases and liquids with largely 
unknown physical properties. Here, we address the important ability of supercritical fluids to sustain 
collective excitations. We directly study propagating modes on the basis of correlation functions 
calculated from extensive molecular dynamics simulations, and find that the supercritical system 
sustains propagating solid-like transverse modes below the Frenkel line but becomes devoid of trans¬ 
verse modes above the line where it supports longitudinal modes only. Important thermodynamic 
implications of this finding are discussed. We directly detect positive sound dispersion (PSD) be¬ 
low the Frenkel line where transverse modes are operative, quantitatively explain its magnitude on 
the basis of transverse and longitudinal velocities. PSD disappears above the Frenkel line which 
therefore demarcates the supercritical phase diagram into two areas where PSD does and does not 
operate. 

PACS numbers: 


Most important thermodynamic and dynamical prop¬ 
erties of an interacting system are related to collective ex¬ 
citations. This has been widely ascertained in solids and, 
to a smaller degree, in liquids [11-0. In contrast, there is 
nearly no discussion and understanding of what kind of 
collective modes propagate in the supercritical state of 
matter, how they evolve with temperature and what are 
the physical consequences. This is related to the lack of 
understanding of the supercritical state itself Si . In 
existing view, no difference can be made between gases 
and liquids above the critical point. At the same time, it 
was long suspected that the supercritical state can show 
elements of gas and liquid behavior depending on pres¬ 
sure and temperature. We have recently addressed this 
point ii and proposed that a distinction between gas¬ 
like and liquid-like states can be formulated in a rigorous 
sense. We have put forward the concept of the Frenkel 
line (FL) demarcating the gas-like and liquid-like behav¬ 
ior at arbitrarily high pressure and temperature above 
the critical point. Crossing the line results in the qual¬ 
itative change of particle dynamics: particle motion is 
purely diffusive above the line as in a gas, and combines 
diffusive and oscillatory motion below the line as in a liq¬ 
uid. Mathematically, the FL can be defined by the disap¬ 
pearance of the minima of the velocity auto-correlation 
function on temperature increase (or pressure decrease), 
corresponding to the gas-like monotonic decay 0. 

At first glance, the above qualitative change of super¬ 
critical system behavior is dynamical only. Yet crossing 
the Frenkel line should also have an important thermo¬ 
dynamic consequence, through the pronounced crossover 
of collective excitations. Indeed, the crossover from 
oscillatory-diffusive to purely diffusive dynamics corre¬ 


sponds to T = td i , where r is liquid relaxation time and 
td is the minimal period of vibration of transverse waves. 
Combining this with the proposal that a liquid sustains 
propagating transverse waves with frequency above y, 
the FL corresponds to the complete loss of two trans¬ 
verse modes at all available frequencies. The crossover at 
the FL corresponds to specific heat Cv = 2fcB, which in¬ 
cludes the contribution of the energy of one longitudinal 
mode, NT, and kinetic terms of two transverse modes of 
each 0, where N is the number of particles. The 
state where c„ = 2 /cb is expected to be accompanied by 
the crossover of c„ to the state with only one longitu¬ 
dinal mode remaiiiing, the prediction supported by the 
numerical results [l0|. Consistent with the above dis¬ 
cussion, the disappearance of the minima of the velocity 
auto-correlation function (VAF) coincides with c„ = 2fcB 
0. This strongly suggests that the dynamical crossover 
at the FL is intimately related to an important change of 
thermodynamical properties of the supercritical system 
related to the crossover of collective excitations. 

In view of the current interest in the supercritical state 
and importance of collective excitations, it is essential 
to follow their evolution directly across the FL and dis¬ 
cuss physical implications. This is the primary goal of 
this paper. For the first time, we directly study collec¬ 
tive excitations in the supercritical state on the basis of 
correlation functions calculated from molecular dynamics 
simulations and discuss important implications for key 
thermodynamic properties such as energy and heat ca¬ 
pacity. We then turn to interesting dynamical phenom¬ 
ena such as PSD. We directly detect PSD below the FL 
where transverse modes operate but not above the FL, 
and provide a quantitative description of the effect. 
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FIG. 1: The location of points studied in the present work 
(stars) on the phase diagram. Full and open symbols corre¬ 
spond to the rigid liquid below the FL and non-rigid gas-like 
fluid above the FL. The Frenkel line obtained from VAFs and 
from cv = 2kB rules are also shown. The inset shows spe¬ 
cific heat Cv along the two isotherms studied. Filled and open 
diamonds correspond to T = 1.71 and T = 2.33, respectively. 

We have performed MD simulations of the commonly 
used Lennard-Jones (LJ) system with parameters e = 
119.8 K and cr = 3.405 A corresponding to Argon. We 
have used 4000 particles in the cubic box with periodic 
boundaries. The system was equilibrated for l.OTO® time 
steps, followed by 0.5-10® production steps. Nose-Hoover 
and constant-energy ensembles were used for equilibra¬ 
tion and production runs, respectively. The time step 
was 0.001 in LJ units. We have simulated several tem¬ 
perature and pressure points in the supercritical state 
below and above the FL. Accordingly, we have simulated 
three densities at T = 2.33: p = 0.549, 0.6655 and 1.071 
(here and below, T and p are in LJ units) and four den¬ 
sities at T = 1.71: p = 0.357, 0.601, 0.845 and 1.012. 
The point (p = 0.6655, T = 2.33) is very close to the 
Frenkel line while the points (p = 0.549, T = 2.33) and 
(p = 1.071, T = 2.33) are above and below the FL, re¬ 
spectively. The location of points under investigation is 
shown in Fig. [TJ We refer to two regimes of liquid dy¬ 
namics above and below the FL as “non-rigid” gas-like 
fluid and “rigid” liquid, reflecting the system ability to 
sustain solid-like transverse modes. 

To directly address the evolution of collective modes, 
we have studied the dispersion of longitudinal and trans¬ 
verse modes and have calculated the correlation func¬ 
tions CL{k,t) = and Crik^t) = 

|A(J 3 ,(k,t) • J,j(-k,t) -b Jy(k,t) • Jy(-k,t)), where 
J(k, t) = the velocity current Q- 

Dispersion curves of longitudinal and transverse exci¬ 
tations are obtained from the location of maxima of 
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FIG. 2: Dispersion curves for T = 2.33 at three densities stud¬ 
ied. The dashed vertical lines mark the boundary of pseudo- 
Brillouin zone. The speed of sound Cs for Ar is from the NIST 
database. 

Fourier transforms and C'r(k,w) respectively. 

The spectra for T = 2.33 at three densities studied are 
shown in Figs. [5] (a)-(c), together with the dispersion 
curves of both longitudinal and transverse excitations. 

We observe that in the non-rigid gas-like fluid at 
p = 0.549, the dispersion curve is close to the Debye 
approximation w = Csfc, where Cg is the adiabatic speed 
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FIG. 3: Crik, oj) at T = 2.33 at three densities studied, show¬ 
ing the peak and propagating transverse excitations at high 
density and their absence at low density. 


FIG. 4: Spectra of transverse excitations at T = 1.71 at four 
densities. Short dashed lines mark the boundary of the first 
pseudo-Brillouin zone. 


of sound. Small upward deviation from the linear depen¬ 
dence (positive sound dispersion, PSD) is observed, but is 
within the simulation error bars (the errors that include 
the calculation of the correlation functions CL{k,t) and 
Cxikjt) and taking their Fourier transforms). We also 
observe that the transverse excitations are absent at this 
density within the error. The same behavior is seen at 
p = 0.6655 very close to the FL. A notable difference in 
behavior appears at p = 1.071 which corresponds to the 
rigid liquid below the FL: we clearly observe PSD in case 
of longitudinal modes and a distinct transverse branch 
for the transverse ones. 

In order to emphasize the qualitative change of the 
spectrum of transverse excitations on crossing the FL, we 
show examples of functions CTikyUj) for all three densi¬ 
ties at approximately the same value of k (Fig. [3|). We 
observe that no peak is found in the current correlation 
function at p = 0.549 and 0.6655, while the clear peak is 
seen at p = 1.071 signifying the propagating transverse 
mode with frequency lu = 4.38 ps“^. 

The same picture of evolution of transverse modes 
emerges from the simulation at the different temperature, 
T = 1.71. According to Fig. HI no transverse excitations 
are found above the FL at p = 0.357 and 0.601. At 
higher densities below the FL, the transverse excitations 
are clearly seen: at p = 0.845 and p = 1.012, they start 
from k > 0.44 and from the origin, respectively. 

Our results of direct calculation of transverse modes 
have two important implications for understanding the 
thermodynamic behavior of the supercritical state. First, 
the theoretical calculation implies that the disappearance 
of transverse waves gives Cy = 2k^ as discussed above. 
Combined with the directly calculated evolution of trans¬ 
verse waves and their disappearance at the FL, this im¬ 
plies that Cy of the supercritical system is 2fcB at the 
FL. This value is consistent with the Cy calculated from 


our MD simulations at the FL, as was the case in the 
previous work Q. Second, the disappearance of trans¬ 
verse modes at the FL implies that the evolution of the 
supercritical system above the FL is governed by the re¬ 
maining longitudinal mode only, with the particle mean 
free path defining the minimal wavelength possible. This 
is a mechanism different from that governing the system 
vibration energy below the FL. Therefore, we expect to 
find the crossover of c„ at the FL between two different 
behaviors. Such a crossover was seen in the previous MD 
simulations [l^ . 


In addition to thermodynamic implications above, the 
disappearance of transverse waves at the FL has an im¬ 
portant consequence for the phenomenon of PSD men¬ 
tioned above. The physical origin of the PSD has re¬ 
mained controversial, including understanding relative 
contributions of the above Frenkel mechanism and other 
effects such as disorder. Experimentally, PSD has been 
well-studied in viscous liquids where the system starts 
sustaining rigidity at MHz frequencies (see, e.g.. Ref. 
[H). Here, PSD is seen at fairly large wavelengths at 
which the liquid can be considered as a homogeneous 
medium. It is generally agreed [l2| that in this range 
of frequencies, PSD originates from the Frenkel mecha¬ 
nism: at high frequency where liquid shear response and 
shear modulus become non-zero, the propagation velocity 
crosses over from its hydrodynamic value v = 


— to the 

p 


solid-like elastic value v = y -^, where B and G are 

bulk and shear moduli, respectively. At smaller wave¬ 
lengths approaching the length of medium and short- 
range order, the wave feels structural inhomogeneities, 
and disorder of liquids and glasses starts to affect the dis¬ 
persion relationship. PSD, with the relative magnitude 
of few percent, was observed in a model harmonic glass 
and attributed to the “instantaneous relaxation” due to 
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fast decay and dissipation of short-wavelength phonons 
in a disordered system 13|. Later work demonstrated 
that starting from mesoscopic wavelengths, the effective 
speed of the longitudinal sound can both increase and de¬ 
crease [IMi. Different mechanisms and contributions 
to PSD were subsequently discussed 01I3- The instan¬ 
taneous relaxation is likely to be significant close to the 
zone boundary [l^, although large PSD in silica glass 
may be related to the effect of mixing with the low-lying 
optic modes. In water, PSD was discussed on the basis 
of coupling between the longitudinal and transverse ex¬ 
citations, and it was found that the onset of transverse 
excitations coincides with the inverse of liquid relaxation 
time [l3, Hi , as predicted by Frenkel. 

Our finding of the presence and absence of transverse 
modes below and above the FL has two important impli¬ 
cations for understanding the phenomenon of PSD. First, 
we see in Fig. [5] that significant PSD does not exist at 
small density and where no transverse modes propagate, 
but emerges in the state where transverse modes kick in 
at higher density. This is exactly the behavior predicted 
by the Frenkel mechanism above. 

Second and perhaps more importantly, we can quan¬ 
titatively address the magnitude of the PSD in the su¬ 
percritical state. A well-known result for solids is that 
longitudinal and transverse velocities can be computed 


as vi = 


B+±G 


and Vt = \/^ correspondingly where 


G is the shear modulus and B is the bulk one. The 
same expressions can be applied to liquids which sustain 
shear stress at high frequency in the solid-like regime 211. 


Combining the expressions for Vh = 
write 


-, Vt and vi, we 


2 2,^2 
D — 


( 1 ) 


Using our MD data, we have taken vi and Vt from the 
dispersion curves at k points where the observed PSD 
is maximal and where uj{k) is in the quasi-linear regime 
before starting to curve at large k. We have used two 
points on the phase diagram: {p,T) = (1.071,2.33) and 
(1.012,1.71). The velocities were taken at fc = 0.595 
for the former and k = 0.47 for the latter point. 
Using the observed Vh, we have calculated vi using Eq. 
(HD. Comparing the actual observed vi with the calcu¬ 
lated one in Table U we find that the Frenkel mechanism 
describes the PSD well. 

We have complemented Table U by the data for vi, Vt 
and from recent experiments for Fe, Cu, Zn and Sn 
2^ 24|. As before, we have used k points where the ob¬ 
served PSD is maximal and where w(/c) is in the quasi- 
linear regime. For Fe, Cu, Zn, we consider the following 
k points: k = 7.7 nm“^, (first point on the transverse 
branch), k = 7.8 nm”’^ (second point on the branch) 
and k = 8 nm“^ (second point on the branch). For Sn, 



Vh 

[m/s] 

Vt 

[m/s] 

vi (obs.) 
[m/s] 

vi (calc.) 
[m/s] 

LJl 

1860 

750±90 

2110T60 

2060±70 

LJ2 

1590 

735±90 

1870±60 

1805±70 

Fe 

3800 

1870T50 

4370T30 

4370±50 

Cu 

3460 

1510±50 

3890±30 

3875±50 

Zn 

2780 

1620±50 

3330±30 

3350±50 

Sn 

2440 

1220T150 

2890T30 

2820±150 


TABLE I: Comparison of predicted and observed vi. The 
data for LJ systems is from our current MD simulations. LJl 
corresponds to (p = 1.071, T = 2.33) and LJ2 to (p = 1.012, 
T = 1.71). The data for Fe, Cu, Zn and Sn are from Refs. 


large PSD is seen at about k = 3.3 nm“^ at the second 
point on the longitudinal branch. To find Vt at this k, we 
extrapolated the higher-lying transverse points to lower 
k while keeping them parallel to the simulation points, 
yielding vt = 1220 m/s. Similar to LJ systems, we ob¬ 
serve perfect quantitative agreement between calculated 
and observed vi in Table HI 


Coupled with the direct evidence of transverse modes 
below the FL, the quantitative agreement between pre¬ 
dicted and observed vi implies that the Frenkel mecha¬ 
nism of PSD quantitatively accounts for both simulated 
and experimental data in the wide range of k spanning 
at least a half of the first pseudo-Brillouin zone. 


Our observation regarding propagating modes is re¬ 
lated to the important point of their detection. Theoreti¬ 
cally, solid-like modes were predicted by simply observing 
that the system behaves like a solid above the frequency 


25l| . Therefore, we have detected propagating waves 


by direct analysis of the current correlation function, the 
procedure similar to that employed in the experiments. 
A notably different result can be found if a scheme ex¬ 
trapolating the hydrodynamic into the solid-like regime 
is employed instead. To compare the direct results with 
one such scheme, the generalized collective modes (CCM) 
approach [26|, we plot the dispersion curves at the same 
temperature and density conditions as in the GCM study 
(see Fig. 7 in Ref. 27|) in Fig. [SJ We observe that the 
GCM extrapolation gives reasonable results for small k 
only but not for larger k: the GCM systematically overes¬ 
timates the excitations frequency in both rigid and non- 
rigid regimes. As a result, the unphysical PSD in the 
non-rigid fluid is observed in the GCM calculation. A 
known disadvantage of the generalized hydrodynamics is 
related to the phenomenological and empirical nature of 
the method [l^, [2^, the general problem that involves 
more specific issues in the case of the GCM scheme. The 
choice of variables in the GCM method is largely ad-hoc. 
Next, the GCM method is internally inconsistent and 
violates the sum rules because the sum rules can only 
be exactly fulfilled if an infinite number of exponentials 
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■ p=0.601 DC 
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FIG. 5: Dispersion curves of LJ fluid (a) above and (b) be¬ 
low the Frenkel line at T = 1.71. Results from direct MD 
calculations and GCM approximation from are shown. 
Short dashed lines crossing the data curves show the location 
of pseudo-Brillouin zone. The straight lines show the Debye 
dispersion law ui = Cak where the adiabatic speed of sound is 
taken from the NIST database [^ . 


is used (see also 0)- Finally, the GCM extrapolation 
incorrectly describes the important region of large k in¬ 
cluding the crossover to the single-particle regime (see, 
for example, eq. 21 in ll|). This is apparent in Figure 
[S] showing failure of the GCM method at large k. 

In summary, we have addressed the propagation of col¬ 
lective excitations in the supercritical state, and directly 
detected their strong crossover. We have found that the 
supercritical system sustains propagating solid-like trans¬ 
verse mode below the Frenkel line but becomes devoid of 
transverse modes above the line where it supports longi¬ 
tudinal mode only. Our results quantitatively explain the 
phenomenon of PSD originating from transverse modes 
in the supercritical state below the Frenkel line. We have 
directly ascertained that the line demarcates the super¬ 
critical phase diagram into two regions where PSD does 
and does not operate. 
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